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Abstract 



The last decade saw the advent of increasingly realistic epidemic models that leverage on the availabil- 
F^ ity of highly detailed census and human mobility data. Data-driven models aim at a granularity down 

^■^ to the level of households or single individuals. However, relatively little systematic work has been done 

1/^ to provide coupled behavior-disease models able to close the feedback loop between behavioral changes 

triggered in the population by an individual's perception of the disease spread and the actual disease 
spread itself. While models lacking this coupling can be extremely successful in mild epidemics, they 
obviously will be of limited use in situations where social disruption or behavioral alterations are induced 
in the population by knowledge of the disease. Here we propose a characterization of a set of proto- 
typical mechanisms for self-initiated social distancing induced by local and non-local prevalence-based 



^ information available to individuals in the population. We characterize the effects of these mechanisms 



in the framework of a compartmental scheme that enlarges the basic SIR model by considering separate 
behavioral classes within the population. The transition of individuals in/out of behavioral classes is cou- 
pled with the spreading of the disease and provides a rich phase space with multiple epidemic peaks and 
tipping points. The class of models presented here can be used in the case of data-driven computational 
approaches to analyze scenarios of social adaptation and behavioral change. 
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C^ Introduction 

^D Understanding human behavior has long been recognized as one of the keys to understanding epidemic 

^"^ spreading [1^ ^ , which has triggered intense research activity aimed at including social complexity in 

. . epidemiological models. Age structure [3j, human mobility [4-16 and very detailed data at the individual 

^ level [TT-'W are now incorporated in most of the realistic models. However, much remains to be done. 

Models based on social mobility and behavior [2H f22 have shown to be valuable tools in the quantitative 

analysis of the unfolding of the recent HlNl pandemic [2TJ [22], but it has become clear that societal 

C^ reactions coupling behavior and disease spreading can have substantial impact on epidemic spreading 

[2^ ^3 thus defining limitations of most current modeling approaches [24 . Societal reactions can be 

grouped into different classes. First, there are changes imposed by authorities through the closure of 

schools, churches, public offices, and bans on public gatherings [25'-'28^. Second, individuals self-initiate 

behavioral changes due to the concern induced by the disease [29-36^. Behavioral changes vary from 

simply avoiding social contact with infected individuals and crowded spaces [37] to reducing travel [i38| [39] 

and preventing children from attending school. In all cases we have a modification of the spreading process 

due to the change of mobility or contact patterns in the population. In general, these behavioral changes 

may have a considerable impact on epidemic progression such as the reduction in epidemic size and delay 

of the epidemic peak. 

Several studies have been carried out in order to evaluate the impact and role that organized public 
health measures have in the midst of real epidemics [26H28] . However, only a few recent attempts have 
considered self-induced behavioral changes individuals adopt during an outbreak in order to reduce the 



risk of infection. In some approaches individual behaviors were modeled by modifying contact rates 
in response to the state of the disease [27l [29l [30l |36l [40] . In others new compartments representing 
individual responses were proposed [STJ [33l [35] . Finally, in some studies the spread of information in the 
presence of the disease was explicitly modeled and coupled with the spreading of the disease itself [32] . 
However, we are still without a formulation of a general behavior-disease model. 

In this study we propose a general framework to model the spread of information concerning the epi- 
demic and the eventual behavioral changes in a single population. The emergent infectious diseases that 
we consider throughout the manuscript resemble the natural history of an acute respiratory infection 
with a short duration of infectiousness and have mild impact on the health status of individuals in that 
healthy status is recovered at the end of the infectious period. We modify the classic susceptible-infected- 
recovered (SIR) model [41] by introducing a class of individuals, 5^, that represents susceptible people 
who self-initiate behavioral changes that lead to a reduction in the transmissibility of the infectious dis- 
ease. In other words, this class models the spread of 'fear' associated with the actual infectious disease 
spread [35[ [42] . Individuals who fear the disease self-initiate social distancing measures that reduce the 
transmissibility of the disease. The spread of fear depends on the source and type of information to which 
individuals are exposed [32l [43] . We classify the general interaction schemes governing the transitions 
of individuals into and out of S^ by considering behavioral changes due to different information spread- 
ing mechanisms, i.e., belief-based versus prevalence-based and local versus global information spreading 
mechanisms. We provide a theoretical and numerical analysis of the various mechanisms involved and un- 
cover a rich phenomenology of the behavior-disease models that includes epidemics with multiple activity 
peaks and transition points. We also show that in the presence of belief-based propagation mechanisms 
the population may acquire a collective 'memory' of the fear of the disease that makes the population 
more resilient to future outbreaks. This abundance of different dynamical behaviors clearly shows the 
importance of the behavior-disease perspective in the study of realistic progressions of infectious diseases 
and provides a chart for future studies and scenario analyses in data-driven epidemic models. 

Methods 
Epidemic model and basic assumptions 

In order to describe the infectious disease progression we use the minimal and prototypical SIR model. 
This model is customarily used to describe the progression of acute infectious diseases such as influenza 
in closed populations where the total number of individuals N in the population is partitioned into the 
compartments 5'(t), I{t) and R{t), denoting the number of susceptible, infected and recovered individuals 
at time t, respectively. By definition it follows A^ = S{t) + I{t) + R{t). The model is described by two 
simple types of transitions represented in Figure ([l]). The first one, denoted by 5 ^ /, is when a 
susceptible individual interacts with an infectious individual and acquires infection with transmission 
rate /3. The second one, denoted by / ^ i?, occurs when an infected individual recovers from the disease 
with rate /i and is henceforth assumed to have permanent immunity to the disease. The SIR model is 
therefore described by the two following reactions and the associated rates: 

5 + / A 2/, (1) 

I ^R. (2) 

While the I ^ R transition is itself a spontaneous process, the transition from S ^ I depends on the 
structure of the population and the contact patterns of individuals. Here we consider the usual homo- 
geneous mixing approximation that assumes that individuals interact randomly among the population. 



According to this assumption the larger the number of infectious individuals among one individual's con- 
tacts the higher the probability of transmission of the infection. This readily translates in the definition of 
the force of infection in terms of a mass action law [44 , Xs^i = (3I{t)/N that expresses the per capita rate 
at which susceptible individuals contract the infection. In order to simulate the SIR model as a stochastic 
process we can consider a simple binomial model of transition for discrete individuals and discrete times. 
Each member of the susceptible compartment at time t has a probability Xs^iAt = PAtI{t)/N during 
the time interval between t and t -\- At to contract the disease and transfer to the infected state at time 
t + At, where At is the unitary time scale considered that we have set to At = 1 day in simulations. As 
we assume to have S{t) independent events occurring with the same probability, the number of newly 
infected individuals /+ generated during the time interval At is a random variable that will follow the 
binomial distribution Bin [S{t), PAtI{t)/N]. Analogously, the number of newly recovered individuals i?+ 
at time t + At is a random variable that will obey the binomial distribution Bin [/(t),/iAt], where the 
number of independent trials is given by the number of infectious individuals I{t) that attempt to recover 
and the probability of recovery in the time interval At is given by the recovery probability /iAt. In this 
processes we recognize that the stochastic variables define a Markov chain [45l |46] of stochastic events 
{5'(t), /(t), R{t) : t = 0, At, 2 At, . . .} in which the current state of the system is determined only by the 
state of the system at the previous time steps. Formally, we can indeed write the following Markov chain 
relations: 

^(t + At) = 5(t)-/+, 

/(t + At) = /(t) + /+-i?+, (3) 

i?(t + At) = R{t)^R^. 

These equations can be readily used to simulate different stochastic realizations of the epidemic events 
with the same basic parameters and initial conditions. This allows us to analyze the model's behavior 
by taking into account statistical fiuctuations and noise in the epidemic process. The equations can also 
be translated into the standard set of continuous deterministic differential equations describing the SIR 
model by using expected values as 



d,S{t) = -/3^(t)^, 

dtl{t) = -^I{t)^l3S{t)^-^, 
dtRit) = III it). 



The crucial parameter in the analysis of single population epidemic outbreaks is the basic reproductive 
number Rq^ which counts the expected number of secondary infected cases generated by a primary infected 
individual. Under the assumption of homogeneous mixing of the population the basic reproductive 
number of the SIR model is given by 

Ro = -. (4) 

By the simple linearization of the above equations for I /N <C 1 it is straightforward to see that in the 
single population case any epidemic will spread to a nonzero fraction of the population only \i Rq > 1. In 
this case the epidemic is able to generate a number of infected individuals larger than those who recover, 
leading to an increase in the overall number of infectious individuals /(t). The previous considerations lead 
to the definition of a crucial epidemiological concept: the epidemic threshold. Indeed, if the transmission 
rate is not large enough to allow a reproductive number larger than one (i.e., /3 > /i), the epidemic 
outbreak will be confined to a tiny portion of the population and will die out in a finite amount of time 
in the thermodynamic limit of A^ ^ oo. 
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FIG. 1: Schematic representation of the two types of transitions that wih be recurrent in the paper. In panel 
(A) we show the first in which individuals in compartment X interact with individuals in class y, represented 
by the small square, becoming Y themselves. In general the compartment inducing the transition of individuals 
in X could be any other compartment in the model, e.g. M, different from the end-point of the transition. We 
assume the homogeneous mixing of the population so that the rate at which an individual in X interacts with 
individuals in Y and changes status is simply given by the product of prevalence Y/N of Y and the transmission 

rate /3, /3Y/N. This type of reaction can be written as X + y — > 2Y. In the case of the SIR model X — S and 
Y — I. In panel (B) we show the second type. This is a spontaneous transition with rate fi in which an individual 
in compartment Y spontaneously moves to compartment Z. These types of reactions can be written as y — > Z. 
In the SIR model Y = I and Z = R. 

In the following we will use binomial stochastic processes to simulate numerically the progression of 
the epidemics and we will use the continuum limit to provide the analytical discussion of the models. 

Coupling epidemic spreading and behavioral changes 

We need to classify the source and type of information concerning the disease that people use to 
conduct their behavior in order to model the coupling between behavioral changes and the disease spread. 
In other words, while the disease spreads in the population, individuals are exposed (by local contacts, 
global mass media news, etc.) to information [32] on the disease that will lead to changes in their 
behavior. This is equivalent to the coupled spread of two competing contagion processes [32l |33l [35] : 
the infectious disease and the 'fear of the disease' contagion processes. The fear of the disease is what 
induces behavioral changes in the population. For this reason we will assume that individuals affected 
by the fear of the disease will be grouped in a specific compartment S^ of susceptible individuals. These 
individuals will not be removed from the population, but they will take actions such as reducing the 
number of potentially infectious contacts, wearing face masks, and other social distancing measures that 
change disease parameters. In the following we will consider that self-induced behavior changes have the 
effect of reducing the transmission rate of the infection, introducing the following reaction: 

S^ + I^ 21, (5) 

with < r/3 < 1 (i.e., r/3/3 < f3). The above process corresponds to a force of infection on the individuals 
affected by the fear contagion A^f^j = rppi{t)/N . The parameter r^ therefore modulates the level 
of self-induced behavioral change that leads to the reduction of the transmission rate. As the scope of 
the awareness of the disease or of the adopted behavioral changes is avoidance of infection, we assume 
that individuals in the S^ compartment relax their behavioral changes upon infection and return back 
to their regular social behavior. While the above modeling scheme is a straightforward way to include 
social distancing in the system, a large number of possible scenarios can be considered in the modeling 
of the contagion process that induce susceptible individuals to adopt self-induced behavioral changes and 
transition to the state S^ . In particular we consider three main mechanisms: 
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FIG. 2: Schematic representation of the third type of interaction discussed. In this case the transition into 
compartment Y is based on the absolute number of the individuals in the compartment (shown by the small 
square). In general the inducing compartment could be different (e.g. M) than the end-point of the transition. 

• Local, prevalence-based spread of the fear of the disease. In this scenario we assume that 
susceptible individuals will adopt behavioral changes only if they interact with infectious individuals. 
This implies that the larger the number of sick and infectious individuals among one individual's 
contacts, the higher the probability for the individual to adopt behavioral changes induced by 
awareness/fear of the disease. The fear contagion process therefore can be modeled as 

S^I ^S^ ^I, (6) 

where in analogy with the disease spread, Pf is the transmission rate of the awareness/fear of the 
disease. This process defines a transition rate for the fear of the disease that can be expressed by 
the usual mass-action law A;^^^ = f3Fl{t)/N. 

• Global, prevalence-based spread of the fear of the disease. In some circumstances, in- 
dividuals adopt self-induced behavioral changes because of information that is publicly available, 
generally through newspapers, television, and the Internet. In this case the local transmission is 
superseded by a global mechanism in which the news of a few infected individuals, even if not in 
contact with the large majority of the population, is able to trigger a widespread reaction in the 
population. In this case the fear contagion process is not well represented by the usual mass action 
law and has to be replaced by 

0F^-§ ^ XiUs^ = Ml - e-''^% (7) 

where < S < 1. Figure ^ shows the schematic representation of this. 

For small values of S we have a pseudo mass action law [44] of the first order in 5: 

Pp(^1 _ e-*^(*)) = Pp [5Iit) + 0(5^)] . (8) 

The above contagion process acts on the whole population even in the case of a very limited number 
of infectious individuals and the parameter 5~^ identifies the characteristic number of infected 
individuals reported by the news above which the fear spreads quickly in the population similarly 
to a panic wave. 

• Local, belief-based spread of the fear of the disease. In addition to the local prevalence-based 
spread of the fear of the disease, in this case we assume that the fear contagion may also occur by 
contacting individuals who have already acquired fear/awareness of the disease. In other words, 
the larger the number of individuals who have fear/ awareness of the disease among one individual's 
contacts, the higher the probability of that individual adopting behavioral changes and moving into 
the S^ class. The fear contagion therefore can also progress according to the following process: 

S + S^^^2S^, (9) 



where the transmission rate is a^F^ with a modulating the ratio between the transmission rate by 
contacting infected individuals and contacting individuals with fear of the disease. The transition 
rate is defined by the mass-action law ^^g^^^F = al3FS^(t)/N. 

The fear/ awareness contagion process is not only defined by the spreading of fear from individual to 
individual, but also by the process defining the transition from the state of fear of the disease back to 
the regular susceptible state in which the individual relaxes the adopted behavioral changes and returns 
to regular social behavior. We can imagine a similar reaction S^ ^ 5 on a very long time scale in 
which individuals lose memory of the disease independent of their interactions with other individuals 
and resume their normal social behavior. This would correspond to spontaneous recovery from fear as 
proposed by Epstein et al. |35j. However, our social behavior is modified by our local interactions with 
other individuals on a much more rapidly acting time-scale. We can therefore consider the following 
processes: 

SF ^S^2S (10) 

and 

S^ ^R^S^R. (11) 

We can then define two mass-action laws: X^f^s = fiFS{t)/N and X^f^^ = fiFR{t)/N. These mimic the 
process in which the interaction between individuals with fear and without fear, susceptible or recovered, 
leads the individual with fear to resume regular social behavior. Both processes, occurring with rate /i^?, 
tell us that the larger the number of individuals who adopt regular social behavior among one individual's 
contacts, the higher the probability for the individual to relax any behavioral changes and resume regular 
social behavior. The two interactions translate into a unique mass action law: A^^.^^^ + X^p^^ = 
XsF^s = l-i^F{S{t) + R{t))/N. The fear contagion process is therefore hampered by the presence of large 
numbers of individuals acting normally. The spreading of fear is the outcome of two opposite forces acting 
on society, but is always initially triggered by the presence of infectious individuals [271 EH 133 133 I3S] • In 
Table [l] we report all the infection and recovery transitions for the disease and fear contagion dynamics 
and the corresponding terms and rates. We will use those terms to characterize different scenarios 
of interplay between the information and disease spreading processes. Unless specified otherwise the 
numerical simulations will be performed by individual-based chain binomial processes [45] |46] in discrete 
time and the analytical discussion will consider the continuous deterministic limit. In the comparison 
between the analytic conclusions with the numerical simulations we will always make sure to discuss 
the differences due to stochastic effects such as the outbreak probability at relatively small values of the 
reproductive number Rq. In the following discussion Rq will refer to the basic reproductive number of 
the SIR model unless specified otherwise. 

Results 
Model I: Local, prevalence-based spread of the fear of the disease 

The first model (Model I) we consider is the coupling of the SIR model with local prevalence-based 
spread of the fear of the disease. The coupled behavior-disease model is described by the following set of 
equations: 
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TABLE I: In this table we show all the transitions and their rates used in the three different models. In the last 
column we write the model in which the transition has been used. For example, the first transition 5* + / — > 2/ 
is related to the disease dynamic and has been used in all three models. 
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FIG. 3: Model I. Schematic representation of Model I. 



dtS{t) = -\s^iS{t)-X's^s,S{t) + \s.^sS^{t), 

dtl{t) = -^lit) + Xs^iSit) + Xs^^iS^it), 
dtR{t) = nl{t). 



(12) 



A schematic representation of the model is provided in Figure (Is]). Considering Table [l| we can write 
down all the terms, 
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in which 



dtXi{t) = for V t and Xi e [S,S^,I,R] , (14) 



meaning that the total number of individuals in the population does not change. In acute diseases, the 
time scale of the spreading is very small with respect to the average lifetime of a person, allowing us to 
ignore birth and death processes and the demographic drift of the population. This is also the time scale 
over which it is more meaningful to consider the effect of the spread of behavioral changes. Diseases with 
a longer time scale may be equally affected by behavioral changes emerging especially as cultural changes 
toward certain social behavior - for instance sexual habits in the presence of a sexually transmitted 
disease with a long latency period - but in this case the demography of the system should be taken into 
account. 

To explain the equations we can simply consider the negative terms. In particular the first term of the 



first equation in Eq. (13) takes into account individuals in the susceptible compartment S who through 
interaction with infected individuals become sick. The second term takes into account individuals in the 
susceptible compartment S who through interaction with infected individuals change their own behavior. 
The first term of the second equation takes into account individuals in compartment S^ who through 
interaction with infected individuals become sick. It is important to remember that the transmission 
rate for people in compartment S^ is reduced by a factor r^^ due to the protection that they gain on 
account of membership in this class. The last term in the second equation takes into account people in 
compartment S^ who through interaction with healthy individuals, 5, and recovered ones, R, normalize 
their behavior and move back to compartment S. The first term in the third equation takes into account 
the spontaneous recovery of sick individuals. 

It is natural to assume that in the beginning of the disease spreading process the population is fully 
susceptible except for the infectious seeds, which means that we can set S^{t = 0) = 0. At this point the 
behavioral response is not active yet. If the disease proceeds to spread much faster than fear contagion, 
then the model reduces to the classic SIR with basic reproductive number Rq = /3//i. In this case the 
initial spread is well described by I{t ^ 0) ^ S^{t ~ 0) ~ R{t ~ 0) ~ 0. The number of individuals in 
the compartment S^ is of the same order of infectious and recovered individuals. From the conservation 
of the number of individuals follows S{t ^ 0) ^ N. Since S is the leading order, all the terms in the 
equations like X{t)Y{t) in which both X and Y are different from S can be considered as second order. 
Using this approximation we can linearize the system and reduce the equations to first-order ordinary 
differential equations that are easy to integrate. In particular for S^ we can write 

dtS^{t) = -^Prlit) - ^rS^it), (15) 

which has the following solution: 

^ ^ fi{Ro - 1) + /iF V J ^ ^ 

For Ro > 1 fear will spread in the population since the condition fi{Ro — 1) > — Mf is always satisfied. 

The growth of the fear contagion is due to the spread of the infection in the population. 

When fear spreads much faster than the disease, /3i? ^ /3, everyone quickly becomes scared and our model 

reduces to an SIR model with a reduced reproductive ratio Rq = rf3J3 / ji = r/^i^o that is dominated by 

the characteristics of the S^ compartment. 

By considering both stochastic simulations of the model and direct integration of the equations, we 

explored numerically the intermediate regime between these two limits, i.e. j3f/ P ~ C^(l)- The spread of 

the fear of infection contagion in this regime does not significantly affect the timing of the disease spread, 

as showed in Figure Q. In this figure the stochastic fluctuations are demonstrated by 50 individual 

realizations and compared with the median profiles obtained by considering 5 x 10^ different stochastic 
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FIG. 4: Model I for /xf = 0.5 day'^ , r^ = 0.5, fi = 0.1 day'^ , N = 10^ and Ro = 2. We show the medians 
of /(t), evaluated using 5 x 10^ stochastic runs for the baseline (SIR model without fear of contagion) and three 
realizations of the model for different values of /3f- In particular in panel (A) we show the baseline SIR model 
with the same disease parameters. In panel (B) we set /3f = 1 day~^ . In panel (C) we set Pf = 2.5 day~^ . 
In panel (D) we set /3i? = 5 day~^ . It is clear how the peak time is the same for all the scenarios and how the 
number of infected individuals at peak is reduced as 13 f increases. 



realizations. The deterministic solution of the equation for /(t), obtained by direct integration of the 
equations, is well inside the 95% reference range of our stochastic simulations as shown in Figure ([s]). In 
this region of the model's phase space fear simply produces a mild reduction in the epidemic size. 

By increasing the value of Pf it is possible to find a region of parameters characterized by multiple 
peaks. In Figure (p| we show 50 stochastic runs and the median profile obtained from 5 x 10^ runs for a 
set of parameters associated with multiple peaks. After the first wave of infection individuals leave the 
compartment S^ and return to the susceptible state in which they are less protected from the disease. 
The second wave manifests if the number of infected individuals at this stage is not too small and if 
there is still a large enough pool of individuals susceptible to the infection. A closer inspection of 
the parameter space by numerical integration of the deterministic equations yields very rich dynamical 
behavior. Figure ^ displays the phase diagram of the model on Rq-Pf plane regarding different number 
of disease activity peaks for a set of model parameters. As r/3 increases, the region in which multiple peaks 
are encountered shifts to smaller values of Rq and larger values of Pf- Fixing r/3, increasing values of Pf 
increase the number of infection peaks while an increase in Rq leads to a decrease in the number of peaks. 
It is interesting to note that adding a simple modification to the basic SIR model leads to scenarios 
with more than one peak. This is important not only from a mathematical point of view (existence of 
states characterized by multiple and unstable stationary points in the function I{t)) but also for practical 
reasons; in historical data from the 1918 pandemic multiple epidemic peaks were observed [26H28] . By 
increasing the value of Pf to larger and larger values, the spread of the fear contagion becomes increasingly 
rapid with respect to the spread of the disease. It is natural to think in this regime that the reproductive 
number of the disease is characterized by the S^ class. We then have two different scenarios: 

1. If Vjsp/ 11 > 1, then the epidemic size is given by that of an SIR model with P -^ Pvjs] 



10 



^x10' 




150 200 

time(t) 



300 



FIG. 5: Model I fixing [3f — 0.25 day~ , rp — 0.5, /i = 0.1 day~ and Rq — 2. We compare the solution of 
the deterministic equations (red solid line) with the 95% reference ranges of our stochastic solutions. Here we 
consider 5 x 10^ runs that produced at least an epidemic size of 0.1% of the population (N = 10^). 
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FIG. 6: Model I Multiple waves of infection. Fixing /i = 0.1 day~^ , i^o = 2, /3f = 3 day~^ , fip — 0.1 day~^ , 
N = 10^ and r/3 = 0.1 we show 100 stochastic runs of the infected profiles and the median evaluated considering 
5 X 10^ runs in which the epidemic size is at least 0.1% of the population. 

2. If VjsP / II < 1, then fear completely stops the spreading of the disease. 

This is confirmed in Figure ([8| in which we plot the proportion of recovered individuals at the end of 
the epidemic, which is evaluated by the integration of the deterministic equations. We consider different 
values of ^f and r/3 and hold fixed the other parameters. It is clear that for very large values of ^f the 
spreading of the disease is characterized by the reproductive number r^^Ro- 

At the end of the disease epidemic the system enters the so-called 'disease-free' stage. This region of 
the phase space is described by 



iS,S^,I,R) = iS,S^,0,Roo). 



(17) 
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FIG. 7: Model I Phase diagram of infection waves on Rq-^f plane. We display the regions of parameter space 
on Ro-/3f plane exhibiting different number of disease activity peaks for three different values of r/3 = 0, 0.15, 0.3, 
where we have fixed /i = 0.1 day~^ , fip = 0.1 day~^ and N = 10^. The phase diagram has been obtained by 
numerical integration of the deterministic equations in Eq. (13). 
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FIG. 8: Model I fixing /j^f = 0.5 day~^ , ji — 0.1 day~^ , N = 10^ and Ro — 2 we evaluate the normalized 
epidemic size Roo/N for different values of Pf and r^ through direct integration of the equations. Once the 
product r/sRo is smaller than unity, then the epidemic size goes to as /3f ^ oo. 



This regime can be easily derived by setting I{t) = in the set of Eqs. (13). The system is then reduced 
to 



dtS{t) = ^FS^{t) 



dtS^it) = -firS^it) 



S{t)^R{t) 
N 
S{t) + R{t) 

N 



(18) 



dtl{t) = 0, 
dtR{t) = 0. 
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From the last equation it is clear that R{t) = constant = Roo^ and the first and second equations are 
equivalent. It is then possible to find the solution for S^ and S by using the conservation of individuals. 
In particular the equation to solve is 



dtS^it) = -^rS^it) 

= -^iFS^it) 

By integrating this equation directly it is easy to show that fear disappears exponentially: 



5(t)+i?ool 


N-s^{ty 


N 



S^{t) 



-HFt 



In the stationary state, for i — > oo, the system reaches the disease- and fear-free equihbrium: 

(S,S^,I,R) = {N - Roo,0,0,Roo). 



(19) 



(20) 



(21) 



There is no possibility of an endemic state of fear. Fear can only be produced by the presence of infected 
people. As soon as the infection dies out, fearful people recover from their fear by interacting with all 
the susceptible and recovered individuals and become susceptible themselves. 



Model II: Global, prevalence-based spread of the fear of the disease 



The second fear-inducing process we consider is the spread of the fear contagion through mass-media 
(Model II). In order to increase ratings mass- media widely advertise the progress of epidemics, causing 
even the people that have never contacted a diseased person to acquire fear of the disease. In this 
formulation, even a very small number of infected people is enough to trigger the spread of the fear 
contagion. To model this we consider a pseudo mass-action law [44] in which the number of infected 
people is not rescaled by the total population. Hence the absolute number of infected individuals drives 
the spread. The transition rate peculiar to this model can be written as X^/^^f = Pf [^ — e""^^^^^] . The 
equations describing the system read as 



dtS{t)= - Xs^iS{t)-X's^s-S{t)-X'sUs-S{t)^Xs^^sS^{t), 

dtS^{t)= - Xs.^iS^{t)^X's^s-Sit)^^sUs-Sit)-^s-^sS^{t), 

dtl{t) = -^I{t)^Xs^iS{t)^Xs^^iS^{t), 

dtR(t) = /i/(t). 



(22) 



A schematic representation of the model is provided in Figure ([9|. Considering Table [l| we can explicitly 
introduce all the terms. 



dtS(t) = 


-pSit)^-^-PFS{t)[l-e-"(''>]+^ 


dtS^it) = 


r,/35^(t)^]J)+/3^5(t)[l e-^^W] 


dtl{t) = 


-^^I(t) + (3S{t)^-^ + r^PS^{t)^-^, 


dtR{t) = 


Hl{t), 



^pS^it) 



S{t) + R{t) 



f^pS^it) 



N 
S{t) + R{t) 

N 



(23) 
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L^-H^ 



FIG. 9: Model II. Schematic representation of Model II. The pseudo mass-action law is represented by the 
dashed line. 



yielding that the population size is fixed, 



^^dtXiit) =0 for V t and Xi e [S,S^,I,R] 



(24) 



As in the previous model, if the infection spreads faster than the fear contagion, then the reproductive 
number is simply Rq = P/ja. In the opposite limit it is easy to understand that the reproductive number 
is Rq = VjsRq. In this latter limit, if r^^i^o < 1, then the global prevalence-based spread of fear suppresses 
the spread of the disease. Moreover, in general we will have a reduction in the epidemic size as a function 
of Vjs. The early time progression of S^ is analogous to that of Model I: 



S^{t) 



SpF 



li{Ro -1) + hf 



(e^ 



(fio-l)t 



-MF* 



(25) 



The analogy is due to the fact that as in the first model the transition to 5*^ is related only to the 
presence of infected individuals. Even in this case the condition /j.{Ro — 1) > —/^f is always satisfied so 
that if Ro > 1, then fear can spread in the population. 

Interestingly, there is a region of the phase space in which this model and Model I are equivalent. In both 
models the transition to fear is related only to the presence of infected individuals. In the first model we 
use a mass-action law while in the second we use a pseudo mass-action law. It is possible to relate one of 
the transmission rates of fear to the other by tuning the parameters. Let us focus our attention on small 
values of S. We can approximate the transition rate by 






Let us consider the first order term only, i.e., X^J^^^f 
sion rates can easily be obtained by imposing 



f3F[SI{t)^0{S% (26) 

PF^I{t)- The relation between the two transmis- 



A 






= X 



which leads to 



P'J 



/3'f 



S^S''' 



1 



(27) 



(28) 



where we define ^^ as the rate in the second model, given (3^ in the first. The above relation guarantees 
the equivalence of the two models at the first order on S. In the small S region in which the approxima- 



tion (26) holds. Model I and II are mathematically indistinguishable for suitable values of the parameters. 



which indicates that even in the phase space of Model II we have multi-peak regions. These regions, of 
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course, coincide with the regions in the first model. 

The disease- free equilibrium of this model does not allow for an endemic state of fear, 

(5'oo, S^, loo.Roo) = {N — Roo, 0, 0, Roo), (29) 

as the transition to fear is induced by the presence of infected individuals only. As soon as the epidemic 
dies out the in- flow to the S^ compartment stops, while the out-flow continues to allow people to recover 
from fear. When the number of infected individuals goes to zero, the media coverage vanishes, as does 
the fear it causes. 

Even in this model the effect of fear results in a reduction of the epidemic size. This reduction is a function 
of 6 and of all of the parameters. As 5 increases the transition into fear becomes faster. Since the people 
in compartment S^ are more protected from the disease, the epidemic size inevitably decreases. While 
keeping the value of 5 fixed, increasing Pf reduces the epidemic size and drives it to its asymptotic value. 
The asymptotic value of Roo as a function of Pf depends on the product rf3J3/ ji. If this product is bigger 



than 1, obtained through direct numerical integration of the equations as shown in Figure (10)- A, then 



the asymptotic value is equal to the epidemic size of an SIR model with (5' = Pr^. If the product is 



smaller than 1, obtained similarly through direct integration of the equation as shown in Figure (lO)-B, 
then the asymptotic value is zero; the rate of the spread of awareness is infinitely faster than the spread 
of the disease. This dynamic can be thought as that of an SIR with a reproductive number smaller than 
1. 



Model III: Local, belief-based spread of the fear of the disease 

In this section we introduce the last model (Model III) in which we also consider self-reinforcing fear 
spread which accounts for the possibility that individuals might enter the compartment S^ simply by 
interacting with people in this compartment: fear generating fear. In this model people could develop fear 
of the infection both by interacting with infected persons and with people already concerned about the 
disease. A new parameter, a > 0, is necessary to distinguish between these two interactions. We assume 
that these processes, different in their nature, have different rates. To differentiate them we consider that 
people who contact infected people are more likely to be scared of the disease than those who interact 
with fearful individuals. For this reason we set < a < 1. 
Let us consider the case of the limit in which no infected individuals are present in the population. The 

S^ compartment can only grow through the interaction S ^ S^ — ^ 25*^. It is possible to show that in 
the early stage this can be thought of as an SIS-like model. Let us consider the case in which there are 
no infected individuals and just one individual in the compartment S^ ^ i.e., 6'^(t = 0) = 1. Considering 
this limit, the set of equations of Model III could be written as 



dtl{t) = 0, 
dtR{t) = 0. 

We assume that in this early stage ah the population is almost fully susceptible S{t ^ 0) ^ N. The 
equation for S^ is then 



dtS^{t) = appS^it) - fipS^it) 



■ /3f : 

a 1 



/j.FS^{t). (30) 
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FIG. 10: Model II Reduction of the epidemic size as a function of /Sf for different values of 6 and r^. We fix 
Rq — 2^ ji — 0.4 day~^, jip — 0.5 day~^ and N — 10^. In panel (A) we assume r/3 = 0.6 for which Tj^Rq > 1. 
Increasing the value of /Sf results in an asymptotic value of the epidemic size other than zero. In panel (B) we 
consider r/3 = 0.4. In this case, instead, r^i?o < 1. By increasing the value of /3f the epidemic size is increasingly 
reduced. This effect is stronger for bigger values of 6. The values are obtained by numerical integration of the 
equations. 

This is the typical early-time term for the 'infected' individuals in an SIS model. The spread of fear 
contagion will start if 



1>0. 



flF 



This allows us to define the reproductive number of fear by 

Pf 



Rf 



flF 



(31) 



(32) 



In isolation, the fear contagion process is analogous to the reproductive number of an SIS or SIR model 
with transmission rate a^F- However, in the general case the spread of the fear of infection is coupled 
with the actual disease spread. The complete set of equations is 



dtS{t)= - 


- \s^iS{t) - \'s^s-S{t) - X's'is-S{t) + Xs^^sS^it), 


dtS^{t)= - 


- Xs-^iS^it) + Xi^s-S{t) + X'/is-S{t) - Xs.^sS^it), 


dtl{t) -- 


= -iil{t) + Xs^iS{t) + XgF^jS^it), 


dtRit) -- 


= ^il{t). 



(33) 
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FIG. 11: Model III. Schematic representation of Model III. 



A schematic representation of the model is provided in Figure (11 ). 
Considering Table |l] we can write all of the terms explicitly, 



dtSit) = 
dtS^it) = 



-pS{t)^M _ p^S{t) 



I{t)+aS^{t) 



-rpfiS^it) 



m 

N 



■ pFS{t) 



N 



fipS^it) 



S{t) + R{t) 



Jit) 



N 

m 



fipS^it) 



N 
Sit) + Rit) 



N 



dJit) = -^I(t)+pSit)'-^+r0pS^it) ^ 



dtR{t) = fil{t). 
Also in this model we assume that the population size is fixed, 

^dtXi{t) =0 for V i and Xi e [S,S^,I,R] . 



(34) 



(35) 



If we consider the case in which the disease spreads faster than the fear of it, then the reproductive ratio 
is Rq = PI ji. In the opposite case the reproductive ratio is governed by the compartment S^ so that 
^0^ = ^/3^o and the epidemic size wih be reduced depending on the value of r^. In this latter case, if 
Tj^Rq < 1, then the protection from infection gained in the compartment S^ causes the disease to fade 
out. Following the same linearization strategy shown in previous sections, the early stage of the S^ 
compartment is given by 



S^{t) 



/?f 



Mi?0-1)-MF(J?F-1) 



,M(fio-l)t 



„|[tF(-RF-l)t 



(36) 



Two different regions in the parameter space are then identified: one in which the rate of increase of 
fear is dominated by its own thought contagion process, jj.f{Rf — 1) > fJ^{Ro ~ l)^ and one in which the 
rate of the local belief-based spread is dominated by the disease, /i(i?o ~ 1) > IJ^f{Rf — !)• In the first 
case the fear spreads independently of the value of Rq^ and the epidemic size will be reduced due to the 
protection that individuals gain in the 5*^ compartment. 

The new interaction, although intuitively simple, significantly complicates the dynamics of the model. 
In particular within several regions of the parameter space we observe two epidemic peaks as demonstrated 



in Figure (12). In this figure we plot the medians for two different values of Rp evaluated considering at 
least 5 X ICrruns in which the epidemic size is at least 0.1% of the population. We also show 50 stochastic 
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FIG. 12: Model III Multiple waves of infection. Fixing fip = 0.5 day~^ , r^ = 0.42, a = 0.05, Ro = 2, 
/i = 0.4 day~^ and N — 10^ we show 100 stochastic runs and the medians evaluated considering 5 x 10^ runs for 
two different values of Rf- In panel (A) Rp — 1.2. In panel (B) Rp — 1.4. 



r,= 0.15 



25.0 




10.0 



FIG. 13: Model III Phase diagram of infection waves on Ro-/3f plane. We display the regions of parameter space 
on Ro-/3f plane exhibiting different number of disease activity peaks for three different values of r/3 = 0, 0.15, 0.3, 
where we have fixed /x = 0.4 day~^ , fip — 0.5 day~^ , a = 0.05 and N = 10^. The phase diagram has been obtained 
by numerical integration of the deterministic equations in Eq. (34). 



runs of the model to explicitly visualize the fluctuation among them. This non-trivial behavior can be 
easily understood. Fear reinforces itself until it severely depletes the reservoir of susceptible individuals, 
causing a decline in new cases. As a result people are lured into a false sense of security and return back 
to their normal behavior (recovery from fear) causing a second epidemic peak that can be even larger 
than the first. Some authors believe that a similar process occurred during the 1918 pandemic, resulting 
in multiple epidemic peaks [26H28] . We show in Figure (13) for a set of model parameters the phase 



diagram of the model on Rq-^f plane regarding different number of disease activity peaks as obtained by 
numerical integration of the deterministic equations. The figure should be considered as illustrative as 
we do not have any analytical expression on the sufficient conditions yielding multiple infection peaks. 



Residual collective memory of the disease and its effect on epidemic resurgence 



At the end of the disease epidemic the system enters the disease-free stage. Setting I{t) = and the 
epidemic size to i?oo the set of differential equations becomes 
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dtS{t)= - appSit) 
dtS^{t)= + appSit) 



N 
N 






N- 


-S' 


'it) 




N 




N- 


- s^{t) 



(37) 



AT 



0, 
0. 



dtl{t) 
dtR{t) 

Conservation of the total number of individuals yields the following differential equation for S^: 



dtS^d) = l^F^^ [{Rf - 1){N - S^{t)) - RfRoo] 



with the solution 



We have defined 7 as 



N 



sF=o{t) 



Nj 



Rf -1 + Be-^^^^ ' 



1 = Rf(1 






1, 



(38) 



(39) 



(40) 



where O is a time-independent variable and is a function of the parameters of the model. Interestingly, 
there are two possible disease-free equilibriums. One in which 



where fear dies along with the disease, and the one given by 



7>0 =^ (Soo^S^Joo^Roo) 



(41) 
(42) 



Ro 



Rf 



-,N- 



RpRoo 
Rf-I 



0,i^oo), 



where fear and behavioral changes persist even after the end of the disease epidemic. The condition 
Rf > 1 is necessary but not sufficient in order to have an endemic state of fear, while Rf < 1 is sufficient 
to avoid an endemic state of fear. Unfortunately, the parameter 7 is an implicit function of the whole 
dynamics through the epidemic size Roq. 

The presence of an endemic state, a societal memory of the disease, and associated fear are quite 
interesting features of the model induced by fear's self-reinforcement. In Model I transition to the 
compartment 5*^ is possible only in the presence of infected individuals. However, in this model fear is 
able to sustain its presence in the population if the effective reproductive number of the local belief-based 
spread is larger than unity even if the disease dies out. Unfortunately, this argument cannot be used to 
fix the range of parameters in the phase space with these properties since any linearization at these stages 
of the compartments is not suitable. The possibility of having an endemic state of fear indicates that an 
event localized in time is capable of permanently modifying society with interesting consequences. In the 
case of a second epidemic, the presence of part of the population already in the compartment 5*^ reduces 
the value of the basic reproduction number. To show this let us consider the differential equation for 
the infected compartment / after the re-introduction of the very same infectious virus (meaning that the 
parameters j3 and ji are equal to those of the first infectious disease): 



dtl(t) 



P 



N 



■r^P 



N 



/i 



m. 



(43) 
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The initial condition of the second disease epidemic could be considered to be the disease- free equilibrium 
of the first epidemic. By using Eq. (42) we can express the rate equation of the infected compartment 



during the early stage of the second disease as 



dtl{t) 



Ro 



Re 



N{Rf - 1) 



r^Ro 1 



RfRq 



N{Rf - 1) 



1 



fil{t). 



(44) 



Let us define di = Roo/N as the proportion of recovered individuals at the end of the first epidemic. In 
the case of the re-introduction of the disease into the population we will have an outbreak only if the 
argument in the parenthesis of the above equation is larger than zero, yielding the following condition 
for the reproductive number R^^ of a second outbreak: 



R. 



II 



^^■ 



di{l -VfsRF) 
Rf-1 



> 1. 



(45) 



It is worth noting that the societal memory of the first outbreak increases the resistence in the population 
against the spread of the second outbreak in a non-trivial way. One might be tempted to conclude that 
the new reproductive number is simply provided by the reproductive number of an SIR model with an 
equivalent proportion of removed individuals (1 — (ii)-, but this is not the case as we have to factor in 
the behavioral changes of individuals in the compartment 5^, obtaining 



Ro <(l-^i) 



^ 
M* 



To prove the last inequality we have to show that 



^/3 



di{l -ti^Rf) 
Rf-1 



<1- di, 



or 



di{l -Tf^Rr) 
Rf-1 



<l-di-r^. 



(46) 



(47) 



(48) 



The expressions on both sides of the above inequality are first-order polynomial functions of r/3. For 
r/3 = 1 they assume the same value —di. It is important to stress that in this limit {r^s = 1) the model 
is indistinguishable from the classical SIR. These two functions can only have one common point which 
occurs at r^ = 1. We will consider only the region in which r^ < 1 as assumed in our model. To prove 
our proposition we have to confront the slopes of the functions and show that 



d\RF 
Rf-1 



< 1. 



(49) 



The polynomial with smaller slope will always be below the other in the relevant region 
can be rewritten as 



< 1. Eq. (49) 



di < 1 



1 
Rf 



(50) 



which is always satisfied, provided our assumption 7 > 0. This is an important result that confirms how 
an endemic state of behavioral change in the population reduces the likelihood and impact of a second 
epidemic outbreak. We note that such a state will inevitably fade out on a long time scale. This can be 
modeled with a spontaneous transition S^ -^ S acting on a time scale longer than the epidemic process 
itself. 
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FIG. 14: Model III Reduction of the epidemic size as a function of Rf and Rq. Fixing r^ — 0.4, ji — 0.4 day~^ , 
Ijlf — 0.5 day~^ , N = 10^, and a = 0.05. The three hues are curves of Roo/N as a function of Rf, keeping 
Ro constant. We select three different values of Ro : 1.5,2.5,3 which correspond to solid black, red, and dashed 
lines, respectively. The value Ro = 2.5 is a special case that leads to Rqv^ = 1. It divides the phase space in two 
different regions. All the values of Ro below are characterized by Ror^ < 1. In this case for large values of Rf the 
model is reduced to an SIR with reproductive number i?o^^ below 1 and the epidemic is halted. Interestingly, this 
behavior starts in an intermediate regime of Rf- There is a critical value Rp of Rf above which (i.e., Rf > Rf) 
the epidemic size is zero. This transition happens with a jump, as shown by the solid black line. All the values 
of Ro above 2.5 are instead characterized hy Ror^ > 1. Also in this case the model is reduced to an SIR with 
reproductive number RoVfs for large values of Rf, but in this case this value is above 1. This results in a epidemic 
size that is always non-zero. In this region of parameters no jumps are present (see the dashed line). The values 
shown in the plot are computed through numerical integration of the equations. 



Discontinuous transition in the epidemic prevalence 



A further interesting characteristic of this model resides in the reduction of the epidemic size as shown 
in Figure (14). In this plot we show Roo/N ^ evaluated through direct integration of the equations, as a 



function of Rp and Rq^ keeping fixed the other parameters. In this case the self-reinforcement mechanism 
creates a more complicated phase space that allows for a jump in the epidemic size as Rf increases above 
a critical value R^ (see the black solid line in Figure (14)). This behavior, typical of the first-order phase 
transitions in cooperative systems, signals a drastic change in the dynamical properties of the behavior- 
disease model. If Rf < 1, then obviously the fear of the disease is not able to affect a large fraction of 
the population and the disease spreads as usual in the population, affecting at the end of its progression 
Rqo individuals. If Rf > 1 we face two different scenarios or two different regions of Rq separated by the 



red solid line in Figure (14): 



In the case that i?o^/3 > 1 (i-e., the dashed line in Figure (14)) the generation of a finite fraction of 
individuals in the S^ compartment is not able to halt the epidemic. The behavioral changes are not 
enough to bring the reproductive number below the epidemic threshold and Roo decreases smoothly 
because of the epidemic progress with a progressively lower effective reproductive number. 

If Rqt^s < 1, (i.e., the black solid line in Figure ( |14| ) the individuals that populate the S^ com- 
partment keep the spread of the epidemic below the threshold. In principle, the state Roo = and 
S^ = N would be possible. In general, the process needs to start with infectious individuals that 
trigger the first transitions S -^ S^ and therefore a small number of R^o individuals are generated. 
However, there will be a Rp at which the growth of the fear contagion process is faster than the 
growth of the epidemic with a small Roo - At this point the fear contagion process is accelerated by 
the growth of individuals in S^ while the epidemic spread is hampered by it. The S^ is quickly 
populated by individuals while the epidemic stops, generating a very small number of Roo- This 
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FIG. 15: Model III S^/N and Roo/N for r^ = 0, i?o = 1-5, /x = 0.4 da2/"\ i^f = 0.5 (ia2/"\ A^ = 10^ and 
a = 0.05. The values are obtained by numerical integration of the equations. 

generates a jump in the amount of individuals that experience the infection as a function of Rf- 
This is clearly illustrated by Figure ( 15 ) where the behavior of both quantities R^o and S^ is plotted 



close to the transition point. The value at which the transition occurs also depends on the other 
parameters of the model including Rq and r/3. 

The extremely rich phase space of this model is important for two reasons: i) we have a strong reduction in 
the cumulative number of infected individuals associated with discontinuous transition; ii) in the case of a 
second epidemic the memory of the system shifts the reproductive number towards smaller values. These 
are very interesting properties of the model due to the self-reinforcing mechanism that clearly creates 
non-trivial behaviors in the dynamics. We have tried different analytical approaches to get more insight 
into the phase transition. Unfortunately, the discontinuous transition is triggered by model behavior out 
of the simple linearized initial state and it is extremely difficult to derive any closed analytical expression. 
An analytic description is beyond the scope of the present classification of behavior-disease models and 
is the object of future work on the model. 



Discussion 



We introduced a general framework with different mechanisms in order to consider the spread of 
awareness of a disease as an additional contagion process. Three mechanisms were proposed. In the 
first, basic model the social distancing effects and behavioral changes are only related to the fraction of 
infected individuals in the population. In the second we modeled the spread of awareness considering 
only the absolute number of infected individuals as might happen in the case that the information the 
individuals rely on is mostly due to mass media reporting about the global situation. Finally, in the third 
model we added the possibility that susceptible people will initiate behavioral changes by interacting with 
individuals who have already adopted a behavioral state dominated by the fear of being infected. This 
apparently simple interaction allows for the self-reinforcement of fear. We have found that these simple 
models exhibit a very interesting and rich spectrum of dynamical behaviors. We have found a range 
of parameters with multiple peaks in the incidence curve and others in which a disease-free equilibrium 
is present where the population acquires a memory of the behavioral changes induced by the epidemic 
outbreak. This memory is contained in a stationary (endemic) prevalence of individuals with self-induced 
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behavioral changes. Finally, a discontinuous transition in the number of infected individuals at the end of 
the epidemic is observed as a function of the transmissibility of fear of the disease contagion. At this stage 
the study of these properties has been mostly phenomenological and we have focused on minimal models 
that do not include demographic changes and spontaneous changes in the behavior of individuals such as 
the fading out of an epidemic over a long time. We should also note that the behavior-disease models 
we have suggested do not take into account the associated costs of social-distancing measures adopted by 
individuals, such as societal disruption and financial burden. A game theoretical approach [29l [30] would 
be well suited in order to account for factors in the decision making process for self-initiated behavioral 
changes. However, more features added to increase the realism of the models inevitably increase their 
complexity. Moreover, the non-trivial dynamic behavior of the models emphasizes the importance of 
calibrating those features by appropriate choices of parameter values. Unfortunately, in many cases we 
lack the data necessary for calibrating the behavioral models. The availability of real- world, quantitative 
data concerning behavioral changes in populations affected by epidemic outbreaks is therefore the major 
roadblock to the integration of behavior-disease models. Any progress in this area certainly has to target 
novel data acquisition techniques and basic experiments aimed at gathering these data. 
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